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SUMMARY 


A  new  two  parameter  family  of  life  length  distributions 
Is  presented  which  is  derived  from  a  model  for  fatigue.  This 
derivation  follows  from  considerations  of  renewal  theory  for 
the  ujnbcr  of  cycles  needed  to  force  a  fatigue  crack  extension 
to  exceed  a  critical  value.  Some  closure  properties  of  this 
family  are  given  and  some  comparisons  made  with  other  families 
such  as  the  lognormal  which  have  been  previously  used  in 
fatigue  studies. 


1.  INTRODUCTION 


It  la  well  known  that  for  the  amount  of  fatigue  data  which  can  usually 
be  obtained  almost  any  two  dimensional  parametric  family  of  distributions 
can  be  made  to  fit  reasonably  well.  In  fact,  in  the  region  of  central 
tendency  the  lognormal,  the  Wcibull,  the  Gamma,  etc.,  can  all  be  fitted  by 
parametric  estimation  and  because  of  the  relatively  small  sample  sizes  hardly  any 
can  be  rejected  by,  say  a  Chi-square  Goodness  of  Fit  test.  However,  when 
it  becomes  a  question  of  predicting  the  "safe  life"  say  the  one  thousandth 
pe.rcentile,  there  is  a  wide  discrepancy  between  these  models. 

For  this  reason  a  family  of  distributions  which  is  obtained  from 
considerations  of  the  basic  characteristics  of  the  fatigue  process  should 
be  more  persuasive  in  its  implications  than  any  ad  hoc  family  chosen  for 
extraneous  reasons.  In  this  paper  we  derive,  using  some  elementary  renewal 
theory,  a  fwo  parameter  family  of  nonnegativo  random  variables  as  an 
idealization  of  the  number  of  cycles  necessary  to  force  a  fatigue  crack 
to  grow  to  a  critical  value.  We  then  examine  some  of  its  relevant 
properties . 
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2.  A  MODEL  FOR  THE  DISTRIBUTION  OF  LIFE 

We  propose  a  basic  framework  and  notation  which  is  similar  to  that 
used  previously  in  (  3  ] .  We  consider  only  standardized  material  specimens 
which  are  subjected  to  fluctuating  stresses  by  a  periodic  loading.  By  a 
load  (or  load  oscillation)  we  mean  a  continuous  unimodal  function  on  the 
unit  interval .  the  value  of  which  at  any  time  gives  the  stress  imposed 
by  the  deflection  of  the  specimen.  Let  be  the  sequence  of 

loads  which  are  to  be  applied  at  each  oscillation  so  that  at  the  1th 
oscillation  load  l ^  is  Imposed.  We  suppose  that  the  loading  is  cyclic 
in  the  sense  that  for  some  m  >  1  and  all 

ljm+i  “  W  £°r  a11  j  *  k  (2‘1'1) 

and  the  loading  is  continuous  so  that  for  alL 

i1+1(°)  -  ^U).  (2.1.2) 

Hence  the  (j+l)8t  cycle  is  the  loading  (tjm)j . fj^|u<). _ 

We  assume  that  fatigue  failure  is  due  to  the  initiation,  growth  and 
ultimate  extension  of  a  dominant  crack.  At  each  oscillation  this  crack  is 
extended  by  some  amount  which  Is  a  random  function  due  to  the  variation  in 
the  material,  the  magnitude  of  the  imposed  stress  and  a  certain  number  of 
the  prior  loads  and  perhaps  the  actual  crack  extensions  caused  by  the 
prior  loads  in  that  eyrie. 

Thus  we  now  make  our  first  assumption 

1*  The  incremental  crack  extension  following  the  application 

of  the  i^  oscillation  is  a  random  variable  with  a  distribution  which 
depends  upon  al!  the  loads  and  actual  crack  extensions  which  have 


preceded  it  in  that,  cycle. 
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The  crack  extension  during  the  (J+1)  cycle  is 


V  J'0'‘ . 


where  the  (possibly  microscopic)  crack  extension  following  the 

load  1^  applied  In  the  i*'  oscillation  of  the  cycle. 

It  follows  from  Assumption  1*.  regardless  of  how  much  dependence 
exists  between  the  successive  random  extensions  per  oscillation  in  each 
cycle,  that  the  random  total  crack  extensions  per  cycle  are 
independent.  Our  notation  will  te 

n 

u  -  V  v,, 

n  ,  j 
j  =  l 

which  has  distribution  function 

Hn(w)  ■  P[U'n  ivl»  for  n»l,2t...  . 

It  follows  from  elementary  probability  considerations,  see  p.  3  69  ,  {6  ] , 
that  the  distribution  of  C,  the  number  of  such  cycles  until  failure,  in 
the  case  failure  is  defined  as  the  crack  length  exceeding  some  fixed 
critical  length  w  for  the  first  time,  is 


PlC^n]  -  1  -  Hn(w). 


(2.2) 


If  we  suppose  that  we  are  dealing  with  a  long  cycle  of  *r.-»  jilat  ions 
which  is  as  complex  as  the  ground-ai r-ground  cycle  in  aeronautical  failure 
studies,  it  might  bo  reasonable  to  assume  *.hat  each  cycle  itself  consist, 
of  a  large  number  of  distinct  phases  of  loading.  Even  though  the  total 
crack  extension  per  cycle  is  the  sum  of  random  variables  which  are  not 
necessarily  Independent  or  identically  distributed,  the  number  of 
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might  be  sufficiently  large  to  make  it  reasonable  to  assume  that 
the  total  crack  extension  per  cycle  is  normally  distributed. 

Thus  we  could  formally  make  a  second  assumption 

2*  The  total  crack  extension  due  to  the  J  cycle  Is  a 

2 

normal  random  variable  with  mean  y  and  variance  o  for  all 

\a  a  possible  eltemsMve  one  could  postulate  that  the  distribution 
of  the  crack  extension  might  be  different  for  the  earlier  cycles  than  it 
would  be  for  the  later  ones.  One  of  the  first  suppositions  might  be  to 
make  the  distribution  of  the  extension  per  cycle  depend  upon  the  size  of 
the  crack  at  the  start  of  the  cycle.  One  such  assumption  and  the  result¬ 
ing  distribution  of  the  total  crack  length  at  the  end  of  n  cycles  is 
presented  in  Section  4,  We  do  not  make  that  assumption  here 
but  instead,  for  reasons  of  simplicity,  proceed  with  the  study  of  the 
implications  of  the  one  above. 
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and 

f.(t)  -  t*4  -  t-**.  (2.7) 

This  two-parameter  family  of  distribution?  is  a  plausible  model  for  the 
distribution  of  fatigue  life.  The  set  of  all  life  distributions  of  the 
form  (2.6)  for  a,8  >  0  will  be  denoted  by  £7  and  most  of  this  paper 
will  be  devoted  to  the  study  of  its  properties.  We  shall  also  refer  to 
the  law  which  has  the  distribution  (2.6)  with  the  somewhat  shorter 
notation  F(a,B) . 

Of  course,  there  are  motivational  derivations  for  other  distributions 
as  well.  The  Weibull  distribution,  which  Is  well  known  for  i  .s  applications 
as  a  life  length  distribution  and  for  fatigue  life  in  particular,  is  obtained 
as  a  special  case  of  the  extreme  value  distributions,  see  p.  302  [7], 

The  Gamma  family  has  also  been  obtained  as  a  distribution  of  life  by 
utilizing  a  model  of  a  bundle  of  strands  which  are  supporting  a  tensile 
load,  see  (4], 

It  is  instructive  to  make  a  comparison  between  the  derivation  of 
the  family  of  distributions  £7  and  an  appropriate  adaptation  of  the 
classical  heuristic  argument,  found  for  example  p.  219,  Cramdr  (5),  as 
it  might  be  used  to  obtain  the  lognormal  distribution  of  the  time  until 
failure  in  fatigue.  Ue  present  this  derivation  in  Section  5. 
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3.  SOME  PROPERTIES  0 T  £7 

Let  T  have  the  distribution  defined  in  (2.6). 

Note  that  T  is  a  two-paramc ter  random  variable  with  6  as  a 

location  pnrametrr  since  it  is  the  median  of  th*  distribution.  (We  show 

later  that  B  is  neither  the  made  nor  the  mean.)  Notice  also  that  a  is 

a  shape  parameter  and  £  a  scale  parameter.  As  we  have  seen  ~  f (T/3) 

is  a  starvdird  normal  random  variable  with  mean  zero  and  unit  variance. 

2 

If  we  let  X  be  ?E(0,  ^*)  we  see  that,  in  distribution, 


2X  -  5(T/B) . 


(3.1) 


If  we  define  the  function  ;■  by 


^(x)  ■  £  1(2x)  for  all  real  x. 


(3.1.5) 


then 


T  -  6<i(X) 


(3.2) 


From  elementary  algebra  ve  find  that 


C* (x )  *  [p  (x)]' 


where 


o  (x)  “  X  +  /x  +1 


(3.3) 


(3.4) 


So  by  (3.2) 


T  -  6(1  +  2X2  +  2x/l+X2  ] 


(3.5) 


where  X  is  *?  E  ( 0 »  ~) .  Hence  we  have  inmeJiatelv 

a 

2 

t(T)  -  6(1  +  y> 


(3.5.1) 
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E(T2)  -  e2 <1  +  2a2  +  ^~) 

var(T)  -  (ae)2(l  +  $f-)  (3.5.2) 

and  we  note  that  for  fixed  a  the  variance  of  T  increases  as  the  scale 
parameter  (median)  8  increases.  This  is  not  true  for  the  lognormal 
distribution  but  empirical  evidence  shows  chat  it  is  true  for  Che  fatigue 
lives  themselves. 


By  noting  that 

4t--x  +  /C7-  p(-x)  (3.5.3) 

P  (  <) 

we  see  that,  whenever  -X  has  the  same  distribution  as  X,  we  have  by 
(3.3),  in  distribution, 

•»«  <3-6> 


Thus  there  follows  immediately  from  (3.2)  the 


Theorem  3.1.  If  T  has  a  fatigue  life  distribution  F(a.B).  In 

£7  then  “  has  a  distribution  in  £7  given  by  F(a,  — ) .  Moreover,  for  any 
i  p 

real  a  >  0,  the  random  variable  aT  has  a  distribution  in  £7  given  by 
F(a,a6). 


It  is  known  that  every  random  variable  with  distribution  defined  by 
(2.2)  for  which  the  are  nonnegative  and  have  densities  which,  are  Fdlya 

frequency  functions  of  order  2,  has  an  increasing  failure  rate,  see  [1]. 
We  now  show  that  the  random  variable  T  does  not  have  this  property  and 
moreover  the  failure  rate  does  not  even  increase  on  the  average  which  is  a 
weaker  cond^cion,  se'  [2].  However,  it  .lust  barely  fails  to  have  this  last 


condition  satisfied. 
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Rentat'k  3.2.  T  does  not  have  a  failure  rate  which  increases  on 
the  average. 

Proof .  Without  loss  of  generality  take  a  «•  B  **  1  and  define 

Qco  -  -mu  -wuctm. 

Writing  Mill's  ratio 

M(0  ■  eU  -  91(01/52'  (f.) 

then  applying  ^'Hospital's  rule  to  Q(t)/t  for  t  -*•  *  and  using  (4.1.1) 
of  the  next  section  we  obtain 

t  '  i-9a  ect)  ]  4  m  2M(e(oi  * 

We  know  that  as  £  ■*■  *  we  have 

M(0  ■  l  -  -y  +  4-+  ocr5). 

C  C 

Hence  .QA--I  m  i  jjUt  Q(l)  •=  A  ^  .  This  proves  our  contention 

that  <joes  not  always  increase  .  j  | 

Actual  numerical  calculation  shows  that  decreases  slowly  for 

t  >  1.64  as  Figure  1  shows.  Since  normal  random  variables  do  have 
densities  which  are  Polya  frequency  functions  of  order  2  the  reason  the 
failure  rate  can  decrease  is  that  our  summands  fail  to  be  nonnegdive  with 
probability  one. 

We  do  not  regard  this  as  being  a  serious  shortcoming  to  our  family  of 

distributions  when  applied  to  fatigue,  any  more  than  the  negativity  of  the 

% 

normal  distribution  vitiates  its  usefulness  when  applied  to  such,  thing.-,  as 
measurements,  which  theoretical!  ,■  cannot  be  negative. 


Figure  1 


Graphs  of  the  density  f,  the  distribution  F,  the  hazard  rate 
and  the  average  hazard  rate  h  where  for  t  >  0 


f<0 


tV*  ,  t 

- —  e*Pt-  J  " 

2t^2v  L 


2t 


1} 


F(t) 


/  f (x)dx,  q(t) 


-lit, L_ 
l-F(t)  • 


h(t) 
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4.  A  SECOND  MODEL  FOR  THE  DISTRIBUTION  OF  LIFE 

In  this  section  we  replace  the  Assumption  2®  by  the  more  general 

one 

g  f 

2'  The  crack  extension  during  the  (n+1)  cycle,  given 

that  the  total  crack  length  was  s  at  the  start  of  the  cycle, 

is  a  normal  random  variable  with  mean  \i  +  fs,  for  some  constant 

2 

6^0,  and  variance  a  for  each  n-1,2,...  , 

By  formula  (2.2),  to  find  the  distribution  of  C  it  is  necessary  and 

sufficient  that  we  find  the  distribution  of  W  ,  called  H  . 

n  n 

By  assumption 

P[Ynfliy  I  VS] 

and  for  n*l 

H.  (y)  *  9J(^~^)  -«  <  y  <  ®.  (4.1) 

i  o 

Now  by  definition,  setting  2  “.x  -_s _ _  _ _ 

Vl(x>  - /"5!(i^>d“ n‘6>-  <4-2> 


We  can  now  prove  the 

Theorem  4, 1  The  total  crack  length  V  ,  at  the  end  of  the 

t  h  2 

n  cycle,  is  normal  with  mean  u  and  variance  a  where 


un  -  £  [ (l+6)n-l J , 


2  2  (l+6)2n-l 

o*o  , -  . 

(1+6) 2-l 


(4.3) 


Proof  by  induction.  The  statement  is  true  for  n=l.  Assume  it  true 


for  it,  then  from  (4.2) 
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n 


where 


x-u  q 

u  “  1+i  v  “  1+4 


But  this  we  recognize  as  the  usual  convolution  of  two  normal  random 
variables  and  hence  the  result  is 


Wx>  - ? 


■'&<) 


Hence  by  simplification  we  find  that 


pn+l  "  "  +  (1+6)V  4+1  "  °2  +  <1+6>2°n 


and  one  checks  that  formulas  given  do  satisfy  the  recursion  relations 
(4.  3).  II 


Strictly  speaking  a  proper  distribution,  analogous  to  that 
obtained  in  (2.6),  which  would  be  the  continuous  extension  of  the 
present  case,  cannot  be  generated  since  P[C  <  «]  <  1. 

For,  one  can  see  that 

u_ -dl 

P[C  <  n]  ) 

n 

and  as  n  -*•  ®  we  have  u  /a  ■*  T~  J (l+6)^-l  < 

n  n  oc 

Of  course  this  would  not  be  of  practical  significance,  since  we  are 

restricting  ourselves  to  situations  where  ^  ^  3,  see  (2.2.1),  hence 

u  _  u 

lim  —  >  3/2/6  for  0  <  6  <  1,  and  lim  ~  >  3  for  any  6. 

G  0 

n  **  n  n 
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5. 


A 


nvtitmtt  rV- 

i.u 


OF 


TltL  LC 


&  L 

In  our  notation  the  i  oscillation  is  caused  by  the  loading  function 
£^,  and  the  crack,  extension  depends  upon  both  the  loading  function  and  the 
size  of  the  crack  already  attained  at  the  time  when  the  oscillation  was 
begun. 


Consider  the  n  loads  which  have  resulted  in  a  crack 

of  size  U  .  then 
n 


w  -  w  ■  ■  s£l— 

n+1  n  6<W  )  * 
n 


assuming  the  incremental  growth  of  the  crack  is  proportional  to  some 

random  functional  Ln+^  of  the  ioadlng  function  imposed  at  that 

time,  and  some  function  of  the  size  of  the  crack  5(U  ).  Now  we  have 

n 


L,  +• • •+  L 


n 

I 

j-1 


<Ywj 


>1)5(Wj_1) . 


Letting  n  become  large  we  then  have  the  left-hand  side  becoming 
asymptotically  a  normal  random  variable,  due  to  the  periodicity  of 
loading,  while  the  right-hand  side  becomes  in  the  limit 


/: 


6(t)dt , 


If  we  take  d(t)  «  1/t  then  Z  becomes  asymptotically  lognormal. 

The  important  distinction  between  this  derivation  and  that  which 
we  previously  considered  is  that  while  a  given  number  of  identically 
distributed  summands  does  approach  normality  as  the  number  gets  large, 
a  random  number  with  large  mean  of  identically  distributed  summands 
does  not  necessarily  do  so.  Moreover,  we  maim. .in  that  what  is  at 
question  here  is  the  number  of  impulses  required  to  exceed  a  certain 
critical  crack  size  anJ  it  is  indeed  a  random  number  of  cvcK-s  which  will 
be  necessary  to  acrrr.pl  ish  this. 


6.  CONCLUSION 


A  derivation,  based  on  plausible  physical  considerations,  for 
a  family  of  distributions  is,  by  itself,  not  a  conclusive  argument 
that  such  a  particular  family  should  be  used  in  life  studies.  No 
family,  however  reasonable  its  derivation,  can  be  accepted  for  use  in 
fatigue  life  studies  until  it  is  confronted  with  actual  fatigue  data 
obtained  under  various  conditions  and  the  distribution  is  shown  to 
represent  adequately  the  life  length  a  which  are  obtained. 

In  order  to  do  this  one  must  have  the  theory  of  estimation 
for  this  family  completed.  The  derivation  of  parametric  estimators 
and  the  ancillary  computing  formulas  for  this  family  will  be  presented 
In  a  later  study.  Also  further  studies  of  the  application  of  this 
distribution  to  the  calculation  of  "safe  life"  will  be  made.  Thus  the 
confrontations  of  this  family  with  actual  data  will  be  carried  out,  to 
provide  the  Justification  for  this  presentation. 
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such  as  the  lognormal  which  have  been  previously  used  in  fatigue 
studies. 
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